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§ 1. Introduction, 

The problem to be considered is the calculation of the polarisation induced 
in a cylindrical rod of finite length by a uniform field of force parallel to its 
axis. 

The only formula previously given appears to be that of Green* for the 
normal component of the polarity at the surface of the rod. In the following 
investigation the method of solution of the general problem of induction by 
series is used, and expressions are given in §§4 and 5 for the mean moment 
of the rod and the induction across the central section. It is assumed that 
the susceptibility of the rod is independent of the applied force ; consequently 
the formulae apply especially to a dielectric rod suspended in a uniform field 
of electric force, this being, in fact, the purpose for which the calculations 
were made. 

In the case ong thin rods in a magnetic field it has been assumed that 
the demagnetising factor might be taken as that of an ellipsoid inscribable 
ill the rod ; but this neglects the fact of the non-uniformity of the field 
within the rod, and various experimental investigations have shown that the 
simple hypothesis is untenable. In § 7 these experimental results are 
discussed in relation to the formulae obtained in the previous sections. 

§ 2. The Series Solution of the General FroUeni, 

Consider a body placed in a magnetic field whose potential is V, and let a 
constant fc denote the magnetic susceptibility of the substance. Let Q be 
the potential of the induced magnetism and ^ the total potential ; then we 
have 

<^- Q + V. 

■^ Green, * Mathematical Papers,' p. 106 ; also Maxwell, * Electricity and Magnetism ^ 
vol. 2, § 439. 
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The components of the intensity of magnetisation I at any point are 
given by 

If dv be the inward-drawn normal at the surface of the body, the problem 
to be solved in finding Q may be expressed in the equation 

Q = J9(I+Q)^-. (1) 



dv ? 



where the integral extends over the surface of the body. 

The exact solution of this problem can only be found in a limited number 
of cases, as, for example, when the surface of tlie body is an ellipsoid. A 
method of approximating to a sohition was given by Beer,* and has been 
further developed by ]Sreumann,f lliecke,+ and Wassmuth.§ The method may 
be regarded physically as one of successive superpositions in the manner 
of Murphy, II Suppose, as a first approximation, that the potential V 
produces an induction Ii equal to — k (d/dx, d/dy, d/dz) V ; then the distribu- 
tion Ii will have a potential Vi, whic^li may in its turn be supposed to induce 
a distribution I2, and so on. The total induction will be given by the 
sum of the series I ; and the success of the method depends upon the 
convergency of this series. 

Defining a series of functions by tho relations 

A 1 = \-rr- —; V2 = --Y^ — ; '^^^^i so on, (2) 

J cv r J vu r 

the first form of series for the induced potential Q, given by Beer, was 

Q = 47r/c Vi -f (47rKyY2 + (^^fcfVi^ + . . . . (3) 

But even if the intermediary potentials decrease in maciuitude, as in most 
practical cases, this seiies for Q is convergent only for a limited I'ange 
of values for k. Other forms of series have been obtained which include 
Beer's series, but give a wider range of convergency ; two such series given 
by Kiecke and ISTeumann are equivalent and are expressed in terms of the 
functions Vi, V2 . . . by 

Q = pVi -r p' (Vi + Vo) + p'^(V, + 2 V2 + V,) + . . . , (4) 

where p = 47r/c/(l+47r/c) = 47r/c/yu,. 

■^ Beer, ' Elektricitiit iind Magiictisiiius,' p. 155. 
+ C. Neumann, ' J3as loo-arithm. Potential,' p. 248. 

I Riecke, ' Wied. Ann.,' vol. 1:5, p. 466, 1881. 

^ Wassnmth, ' Wied. Ann.,' voJ. 51, p. 367, 1894. 

II Murphy, * Electricity,' p. 93, 1833. 
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This series is certainly convergent for all practical cases in which the 
intermediary potentials form a decreasing series, for instance, in the case 
of an ellipsoid placed in a uniform field parallel to an axis ;* for then we 
find 

V= -Ho:; Vi = ^Hx- V2= -/J^)'hx; ; 



where L is the usual constant associated with an ellipsoid and is less 
than 4 tt. 

Forming the functions for series (4) we obtain 



K 






P^ ( Vi -f V2) = ^^-'^- ^^ '^ - LH.^ ; 



p^ ( Vi + 2 V2 + V3) = (4^-L)'^' a LHa^ etc. 
Hence we have 

/^ L /X fJL" J 

But this series is convergent for all values of Ac,and in fact gives the known 
exact solution 

I + kL 

We shall use, then, the series (4) in the following work, in which we consider 
tlie case of a long rod in a uniform field. 

§ 3. Series for a Finite Rod, 

We shall consider the case of a cylindrical rod of length 2 1, with fiat 
circular ends of radius a ; we write Ija = m. 

Suppose the rod placed in a uniform field H parallel to its axis, which 
we take as the axis of x ; then we have 

V=~H^^ (5) 

Our object is to calculate the functions Vi, V2, ... defined in (2). 




* Cy. lliecke, loc. cit, p. 481. 
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Using cylindrical co-ordinates {x, p, 6), we have in general 

Also the element of surface da is equal to p'ds d6\ where ds is equal to dx' 
or dp according as the point {x' y p) is on the curved side or a flat end of the 
rod ; and putting 

R2 = p2^^'2_2/>/c0S(9', 

we have 

47rVn+i {x, p) = d0'[ p' -^ ds e^^(^ '-^)Jo (\E) d\ 

= 27r f />' -J^ ^5 p±M^/-a')J, (Xp) Jo {Xp') dX, (7) 

where the integral with respect to ds extends over the path ABCD, and the 
upper or lower sign is taken in the exponential according as x'—x is 
negative or positive. 
Thus we find 

Vn+l (X, />) = i r(^5^)^, /y p-"^"^^^Jo (V) Jo (VO ^^ 

^i rff r)^, ^/^^/>' Ce-^i^-^^Jo (Xp) Jo (Xp') ^X 

- i f" f??^ ^^^^^' f " ^^ ^"'""^^ Jo O^p) Jo (>^^) ^^ 

]_\dp/p' = a Jo 

-i(Y??^ a^Z^r' re-M-'-)Jo(\/))Jo(\6i)6^\. (8) 

]^\dp/p' = a Jo 

From symmetry in the case under consideration, this reduces to 

Vn+i (x, p) = - r(^] ^ _p'dp' Ce-^' sinh (A^) Jo (Xp) Jo (VO ^^ 

__ f^^?^\ ac^^' f "e-^^' sinh (X^) Jo (Xp) Jo (Xa) dX 

h\^P /p' = « Jo 

- rftV) ^^^^' [" ^"'''* ^osh (X^O Jo (Xp) Jo (Xa) dX, (9) 
Jo\8p A' = « Jo 

Thus at any stage in the evaluation of the successive potentials, V^+i is 
expressed as the potential of a surface distribution due to V^ ; the first 
integral in (9) gives the part due to the flat ends of the rod and the other 
two integrals that due to the curved sides. 

The potential Vi can be deduced from (9), or obtained directly as the 
potential of two uniformly charged circular discs situated at the ends of the 
rod. 
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Using the relation 

fa 
P Jo {^p) dp = aJi {\a)l\ 




we obtain 

.00 



Vi = Ha e-^^ sinh {\x) Jo {Xp) Ji (X^^) dkjX. (10) 

Jo 

§ 4. Oalculation of the Total Moment of the Bod, 

The magnetisation at any point of the rod is equal to — /c (3/3^, djdy, d/dz) <^, 
where <^ is the total potential, and the moment of the rod can be calculated 
by considering the normal component of the polarisation over the surface 
of the rod. Writing M for the total moment, we have 

M = -2/ r27rKp(^i) _rfp-2 C2irKCWo(^) dx, (11) 

where the first integral gives the part due to the ends of the rod, and the 
second integral that due to the sides. 

Now <j> is given in terms of V and the intermediary potentials Yn according 
to the series given in (3) and (4) ; thus from (11) we find that it is required 
to evaluate expressions of the form 

U^ = 21 f 27rp f^™^^ dp-\^2C27rax(^) dx, (12) 

Jo \dx/x=^i Jo \dp/p:=.a 

We shall consider these expressions as far as Ug, and obtain the important 
terms in their expansions in powers of 1/m, where m is the ratio of the 
length of the rod to its diameter, and is assumed to be large. 

Then since 

V = -Ux, 
we have 

U = -2waHfc'H. (13) 

Further, from the value of Vi given in (10) we have 

/»oo 

Ui = ^wa^m I e-^^X-1 cosh (Xl) Ji^ (Xa) dX 

Jo 

• 00 

e~^^ sinh (Xx) Ji^ (Xa) dX, (14) 



Awa^H 1 xdx 

Jo 







If in the second part we perform the integration with respect to x, and put 
fjb for Xa, we find that (14) reduces to 

Ui = 27ra3H C ^^^l^^^i//,-27ra3H f e-^^^^f^dfi, (16) 

Jo /^ Jo /^ 



l^ow we have 

/^^ "^ ~~ ^'^ 



\ 



^JMa^ - A 



D 2 
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Also, since m is large, the important terms of the second integral are 
obtained by substituting the value of the Bessel function for small values of 
the argument ; thus as far as the term m~^ we have 



f 



^-2mf.llJS dll =-i , (16) 



Hence to the same order we have 

U. = 2„'H(±-i) 

\37rm 8mV ' 

Now we have to evaluate U2 in a similar manner ; we have from the- 
general expression given in (9), 

( -~j ) adx' \e-^^ cosh ^\x') J. (Xp) Jo (Xa) dX 

q\CP /a }q 



,i./^\^ifY3li\.v.' 



W^) ccdx' { Xe-^^ cosh (Xx') Jo (Xp) Jo (Xa) dX. (18>> 



And, similarly, from (9) we have 



dp / =a 



^"j + {"(—}] p'dp' f "xe-^' sinh (\x) Ji (Ka) Jo (V) dX. 

(19)> 

It may be noted that if the intermediary potentials V^ be regarded as due- 
to successive surface distributions <7^, then the equations (18) and (19) are 
examples of general relations between consecutive surface densities. 

Using the values of (18) and (19), we have 

U2 = - ^ r2irp (^] dp - 1 r 27rpdp {''(^]p'dp' Cxe'^^' Jo (Xp) eJo (Xp') dX^ 
4- 21 2irpdp ( ^— ; j adx' Xe"^^ cosh (X^') Jo (Xp) Jo {Xa) dX 

— 27raxl'-^] dx + 2 27raxdx\ l-rr-j-) P^Pi Xe~^^ sinh (Xa:?) Ji (X^j) 
Jo \dp /a Jo };,dx li"^ ^ ]^ 

Jo {Xp) dX. (20> 
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Using the formulse 

d\ \p^ (p) Jn (V) Jn (^^) dpz=z(j^ (r), q < r < py 

i io 



Jo 



Jo (Xp) Ji (\d) d\ = -y p < a, 

^ 



we find that (20) reduces to 

U2= ~ f 27raa? (^) &?-4to f Y|^^) p'do' f "e-^^ sinh (XI) J. (Xp') Ji (\(^) — 
Jo \op/^ JoV^A Jo ^ 

+ 47raZ (" /^M adx' j" g-^^ cosh (Xo?') Jo (\a) Ji (Xa) (?X. (21) 

We substitute now for Vi from (10) and write 

X = a^; p z=z a7)\ I = am ; fjL = Xa; /i-i = X^a, 
Hence we obtain 

Jo Jo 

— 47rHa^ 7/^?; e~*^/^i cosh (m/^i) Jo (//.?;) Ji (ytti) c?/ai 

Jo Jo 



JO 

e'-mf, sinh (m/i) Jo (fiv) Ji (ft) ^ 
f^ 



f 



CO 7 /lOO 



27rHa3m ] e-'"'^^ cosh (my.i) Ji^ (//^i) 2£i _ 4:7r'Ra^ j ^-^^'^ sinh (mfi) Ji^ (//,) ^ 



^1 e-^'^isinh (//,i|^) Ji^ (^i) dfn 

Jo 

/»eo 

j 6-*^/^ cosh (yU,^) Jo (/x) Ji (m) ^//'. (22) 

Jo 

The first two terms give the following ; 

^^ -47rHasf " e-^f' sinh (mfi) Ji^ (/^) ^ 

Jo ' ' \ " ^' //'I Jo /^' 

-27rH«^3 [\a7} fV^^^/^jJo (^i7;)Ji(^i)^/ii re-""*'^sinh(m/i)Jo(/^7/)Ji(/^)^, 
Jo Jo Jo f^ 

Also, in the third term of (22) we perform the integration with respect 
to ^ ; then in all the terms of Us we apply the previous method of approxi- 
mation, that is, we consider m large and substitute the lower terms of the 
expansions of the Bessel functions for small values of the argument. Using 
also the relation 

Jo /^ 

and keeping to the order in which the parts of U2 are given above, the 
approximation indicated gives, to the order 1/m^, 
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Us = TrHa^; +— — -- — j^ ^-. + 



= ^Ha2/(__^ +11,). (23) 

Comparing XJ2 with Ui given in (17), we see that they both begin with a 
term of the same order 1/m ; but we have 

Ua + U2=~i-7rHa2^-i. (24) 

24 rii'^ 

Eeturning noM^ to the series (3) and (4) for the induced potential Q, we 
see that in calculating the total moment of the rod in this manner it is of 
advantage to use the series 

Q = ^Vi + tey(Vi-i-V2) + ..., /. = l + 47r^. 

Then, if we denote by 1^. the mean moment per unit volume obtained by 
calculating the moment of the whole bar, we have from (13), (17), and (24) — 
as far as the order indicated — 



Im = kH < 



. _47nc (Jt J_\ MttacV 1 "X .2-\ 

fi \37rm 8mV \ fjL /48m^J • ^ "^^ 



§ 5. The Mean Induction through the Central Section. 

Another quantity which we require to calculate is the mean polarisation 
over the central transverse section of the bar ; if this be denoted by lc> 
we have 

I,= ~-^fV^(^) dp, (26) 

where </> = — Ha;+Q. 

To evaluate the contributions made by Vi and Vg, we have from (10), 

1- [^irp (^ ^/) = ?H [%dp Ce-^^ Jo {\p) Ji {\a) d\ 

m^]^ "^Xdx h "^ a Jo Jo 



ira ,Q 



=:2Il["6^-"''^Ji2(^)^/^/^ 

Jo 



to the order indicated, using previous method of approximation. 
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From the general expression in (9) we have V2 given in the form 

Vg = —Ha [^pdp [ e-^i^ cosh (Xil) Jo (\ip) Ji (Xicc) dXi 

Jo Jo 

e-^^ sinh (\x) Jo (\p) Jq (^p) d\ 

Jo 

4- Ha^ i^dx'i e-^i^ sinh (\ix') Ji^ (\i«.) ^ZXi 

J a: Jo 

1 e~^' 8mh{\x)Jo(Xp) Jo{\a) d\ 

Jo 



-l-Ha^ 



r^aj' 6-^1^ sinh (Xi^x;') Ji^ (Xi^) c^Xi 
Jo 



Hence we have 



/•oo 

j e-^ cosh (\x') Jo i\p) Jo (^a)d\. (28) 
Jo 



3a; Jx = 



-^] =-H« 







/)'^/)' g"-^!^ cosh (\i/) Jo (Mp) Ji (A.!^) ^A-i 
Jo 



[ \e-^^Jo(Xp) Jo {Xp)dX 

Jo 

4- Ha2 c^i^' f e-^i^ sinh (Xio^') Ji^ (Xia) ^\i 

Jo Jo 



\e~^^' Jo (Xp) Jo (X«) ^A'. 



(29) 







The mean value of this taken over the circle x =^ 0, p = a, is given by 

TT^ Jo \CX Jo 

= -2H {'' pdp 6-^1^ cosh (XiO Jo (VO Ji (^1^) ^^1 f e-^^Ji{\a)Jo {\p')dX 

J Jo Jo 



-|-2Ha 



dx < 

Jo Jo 



e-^)^ sinh (\i^'') Ji^ {\ia) dX] 



00 
e-^^' Jo (Xa) Ji (Xa) dX, (30) 







We proceed to evaluate these integrals for m large by substituting the 
expansions of the Bessel functions as before ; then the first integral in (30) 
becomes, to the required order, 

The second integral in (30) is equivalent to 



"^ ' ^ i ^ + Sm^ 1 6 (^4^2 + g j ^^^4 J \ 2m^ ^ \4:a^ "^ 8/ m^j ^^^ 



H 



Jo Jo 1/^1- 






Ji^ (yu-i) Jo (/^) Ji (/^) ^/^ ^/^i- 
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Substituting values of the Bessel functions of fju or fii, or of both, according 
to the exponential term, this reduces to 

Adding this to the part given in (31), we see that the mean value of 
3V2/3^ over the central section begins with a term —23H/32m*, while the 
mean value of 9Vi/9^ begins with a term of the order 1/m^. 

Then, substituting in (26) the values calculated above, we have, as far as 
the order shown, 



L \2m2 4mv ^ ^ 32m*. 

= /scH^ 1 --\-7rK\^-\--~'irK\--~'. >, (33) 

L m"* \ 2 /m*J 



1 6. Application to Dielectric Polarisation. 

If an ellipsoid of specific inductive capacity K is placed in a uniform 
electric field, it is known that the field within the ellipsoid is uniform, and 
the polarisation at any point is given by a formula 

I = /^(H-NI); /c = (K-l)/47r. (34) 

The factor N is a numerical factor independent of I, and can be calculated 
theoretically; for a prolate spheroid, of which m is the ratio of the longer 
to the shorter axis, we have 

N= 4^/--;p^log,(m+^/;^ (35) 

the axis of revolution being along the direction of the field. Also, for an 
infinite rod of circular section, N is zero for a uniform field parallel to the 
axis of the cylinder, and is equal to 2 7r for a field at right angles to 
the axis. 

In the case of a rod of finite length placed in a uniform field, there is no 
longer uniformity of the field within the rod, but as we suppose the specific 
inductive capacity to be independent of the field, we have still an equation 
of the form (34). For a field parallel to the axis of the rod, we have from 
(25) to the first order of 1/m, the ratio of diameter to length, 

N = 16/3Km. (36) 

For a field at right angles to the axis we may still take N" to be equal to 
2 7r. Now, in determinations of the specific inductive capacity of a dielectric 
in the form of a cylindrical rod, it is necessary to measure the couple 
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required to keep the rod at a given angle <^ to the direction of the external 
field. If we put Y for the volume of the rod and use the values of N 
^obtained above, we have 

Couple on rod = kW , ' o^Ti^^^ffT^ir x ^' ^^^ ^ ^^« ^' (^7) 
where k = (K— l)/47r. 



§ 7. The Demagnetising Factor of a Finite Bod, 

When the problem is one of magnetic induction, it is clear that the 
previous results only supply an approximation to a solution ; for, since the 
field is not uniform within the rod, it follows that the magnetic permeability 
jjb varies along the rod, thus complicating the problem to be solved. In fact, 
the demagnetising force will not be expressible in the form NI, where N is 
-a numerical factor independent of L 

On the supposition that there is such a constant IsT, at least approximately, 
for long thin rods, Du Bois* made a comparison between the values of K 
for ovoids and cylindrical rods of the same dimension-ratio m ; the factors 
for the rods were calculated from experiments by Ewing for the longer rods 
.and by Tanakadate for the shorter ones. 

It was pointed out by Mannf that in the former experiments ballistic 
methods were used, but in the latter magnetometric methods. Consequently 
the quantities measured were those denoted by I^ and Ic above ; and the 
factor K" deduced by the two methods will, in general, be different, and may 
be denoted by IST^^ and N"^, being the magnetometric and ballistic values of 
the demagnetisation factor. From the series given in (25) and (33) we have, 
for large values of m, the approximate values 

If„ = J^; N* = ^. (38) 

BenedicksJ further investigated the values of N"^ and ISb for rods and, 
in particular, their variation with the intensity of magnetisation ; the value 
of the demagnetising force IST^I at the centre of the rod was found to increase 
with the field, and then after a certain point to decrease rapidly. 

Using the ballistic method, Lamb§ determined the induction through 
sections taken at intervals along a bar, and was thus able to trace curves 

■^ Du Bois, ' Annalen der Physik,' vol. 46, p. 495, 1892. 
t C. E. Mann, ' Physical Eeview,' vol. 3, p. 359, 1895. 
X Benedicks, * Annalen der Physik,' vol. 6, p. 726, 1901. 
§ C. G. Lamb, 'Phil. Mag.,' vol. 48, p. 262, 1899. 
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showing the distribution of induction and force along the bar and the 
consequent effective distribution of magnetic permeability. 

Searle and Bedford* investigated the demagnetising force /i- at the centre 
of a long cylindrical rod. Determining the values of h for increasing values 
of the external force H, they found that A increases at first more rapidly 
than 27rlc/m^y but after passing a maximum it descends to a minimum 
somewhat less than that value ; the latter result is interpreted as being due 
to hysteresis effects. 

A quantity measured in several experiments was the ratio Im/Ic, giving 
the effective length of the bar. 

Using expressions (25) and (33) above, we have approximately 



With a rod for which m was 25, Benedicks found this ratio increased to 
0*84 at the highest value of H used in the experiment. 

Eeviewing the experimental results, it appears that in the case of magnetic 
induction the facts of the non-uniformity of the force within the rod and 
the consequent variation of effective susceptibility are too important to be 
neglected ; hence, in this case, the calculations of the previous sections 
supply limiting values towards which the quantities approach as the 
magnetisation tends to saturation. 



* Searle and Bedford, ' Phil. Trans.,' A, vol. 198, p. 98, 1902. 



